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GIISAVG is a higk energy laser propagation computer 
prograx fcr grourd-tc-space apelications. Written 2 y el ae 
Ee Hcc¢ge from the Air Force kéeapons Lapcratory, Kirtlane 
AFB, it is cne in a tamily cf fFropagation ccdes addressing 
this application. Specifically, GUISAVG was designed tc 
compute irradiance at the target given a ucdel atmosphere, 
lasé€r device parameters, and simple target enyagement geor- 
¢try. The transmitter induced effects of beam guality and 
jitter are considerec as are tke linear atmosjneric effects 
cE scattering, aiSCEPplacnu, “aid CUED wy -nec wee eer 
kleocaing medel is alsc included. Adaptive optics Ccnpencae 
tion can re applied tith consideration given to iscplanatic 
effects. 
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LT. ENTECLUCTIION 


GCISAVG 1s a simjlified laser propagation progran- Lome 
interded specifically for near vertical ground-to-space 
appiicaticaos utilizirg a £iwed earth—pased  transiatre. 
dicected at a sirgle taryet satellite. The primary furrfose 
ct thé program is te y;rovide irradiance and fluence cn 
tacyet alceng with related propagation data. 

Tre p~regram was written cry Charles B. Hogge from the Air 
Force Weapcrs Lakoratery, Kintland Ait force. acc, 
Alpbustergue, New Mexicc. GUISAVG is one of a family ci 
ground-tc-Szface rropagation pregrams. Other versiors include 
GUISFE {fectprint) ard GUISMIF. GUTSFP computes an engage. 
Jent envelcje based cr user supplied irradiance thresncid 
levels. The propagaticn calculaticn methods are identical to 
GUTSAVG. GUISMIF is a full wave-optics frogram using fast 
Fourier transforms ir the beam fropagaticn computaticns. 

The kasic apryroach used in GUTISAVG is to utilize nucdula- 
tion transfer functicrs to characterize eriects such as Lean 
jguality, jitter, turkulence and to apply these effects at 
the aperture as a sircle phase screen. Tne same aprrcach is 
used in ESF-IV [Ref. 1]. The-general modulation transfer 
functicn, (MIF), usec in this way is described in fcllcwirg 
secticns alcng with the development of each beam degrading 
mechanism. A linearized model is used for thermal blcc#ing 
which is also effectively applied at tne aperture as a fghase 
screer. In the case cf bloo@ing, however, the phase variance 
due tc biconzing is ccaputed, and the Strehl relation is used 
to deteraine the relative irradiance reduction. The effect 
cf thermal rlooming and the other spreading effects are 
combired using beth BSS and multiplicative methods. Finally, 


the average of these two methcds is used to determine the 


10 


Pores yctei trrddiarce reducticn from tne diftracticn 


limited caceé. 
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Ae TEE ETF APPRCACH 


GUISAVG uses MTFs tec apply and tnoen to consolidate V2 
effects cf Fitter, beas quality and curputence: 
zssentiaily, these effects are repiaced by a f~hase screen at 
the aperture which multiplies the initiai complex aperttre 
districruticn by a rald¢m phase distorticn Lactor. 

Lutcuicski and Ytra [Ref. z2] wave developed an Cxpiecee 
sion fcr the average intensity at a point P in the far-field 


jue tc perturbations at the aperture. ine intensity at P is 


2 - 2 
IGS) = (5) fo, (eZ) exp - (#)50] dp 





Za0p 
+ foceme’ (2-1) ex | i ae 
where 
= ees (Zz 
and 
r = 27 (2 


¢ is the vector from the z axis of symmetry in the far-field 
to the ~cint P. U(r) is the aperture distcibution ard By (P 
+Z) is the MIF fer tlre disturbance. ££, and ©, ane txo Gee, 
trary feints in the aperture where the phase perturpaticn is 


measured. P is the distance fetween the points. nom 


PIemeoiciceteasnere, 15 4 Ccnstant. If the phase distur- 
Lance if a KLandonz Gat£sian yaeriakie wath a kuncwn correlation 


Mmiiecich, tren cle NGF Can Fe €xpressed as 


M (8) =<expli(r, —en|> (ea 


Mmieercris Of the structure £Luncticn, 


My(d) = exp [(-4)D,(ri-ra)] (2-5) 


Meenscecccnd integral amequatior 2.1 Lepresents the urrerial- 
Mzemmaperttire MIF. N¢ecriaiizing this term witn the fowrer it 
the arperttre, (Po), tc cause the MTF to be unity at the 

Peg in results in 


2 


sles ie 2 
Co) (#)e.fu, om, Oe kop] 47 (2.0) 


Notirc tke symmetry cf the intensity for a given p and 
expressifrg Hf) as the combined effects ot jitter, rear 


Guality, and turktulence, equation 2.€ can Le rewritten ée@s 
rp) = (%\p,fu. am, ym, (DM. (3) 3 MEE) ap (2 +7) 
Ez a) ae t b a ee eae 


A Fourier-Eessel trarsform has teen used.Note that the FIFs 
ci jitter (hs), ceam guality(M,), and turbulence (M,) have 
Feen sukstituted for Eo: 

ike ESE-IV manual { Ref. 3) contains the development 
above in mere detail. MIFs for the specific effects are 


dGescribrec ir fcllowizc secticns. 


amd 
day 


Be, JPERPAL SEL CO cana 


Tke thermal tlconing nodei in GUTSAVG is tased or tke 
follcsirg jJinearized density perturbation eguation. 


xt 
a (2) 
Vo 


Is 





I ( ae rd t 
x ye rere ae 


=O 
xX exp - flatarto(al] ar] 


OD 
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yo 
“y_ 


TO le 


Eere, a(%) is the atuacspheric absorjtion coerficient at a 
distarce 2 along the Leam, Vy, iS a constant transverse sind 
velocity, ard p, 1s the ambient pressure. ihe exponeneiam 
term represents the tctal extinction due to scatterirg anc 
absoryticn. x'is a ccerstant ci integration. It Lf aeeumen 
that the beam is prcfagated in the positive z directicn a 
distarce 2 and that the wind vector 1S in the positive x 
directicr. The irtensity integral represents tne heating cr 
the atucsphere as it transits the beam [Ref. 4]. Scme 
assumpticns embodied in the akcve eguation are that V,<< Cc, 
. the ilccal sonic velccity, sc that the process represerted 
cccurs at ccnstant pressure, and that the kinetics cf 
absorpticn and ccnversicn tc heat are extremely fast 
[Ref. 5]. | 

Ey applying the Cladstone-Llale relation,the density 
relaticr car be expressed as a Change in the refractive 


index. 


he = (OL (2255 
O09 
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Eguaticn 2.€ can be rewritten as 
x! 
An = eee iat bf roe syiax 
0 
. 00 (22610) 


x exp - fearon] a] 


The chance in the wavefront phase of a beam due to the 
refractive index charge wnern the beam 1S propagated a 
distance eg, is given ry 
oi 4 : 
A¢d = ~- And& (ea) 


r 
0 


m@estatutinsg Eecuaticr 2.10 iantc Z.11 and Changing the limits 


MmmimtcCLadatlon te reflect the Qround-to-Space propa.aticn 


Fath results in x! 
Zig YL 
pean «ACEP ree ae 
o y) 2 oF YJ. ay) ax 
h 
"atm ee) 
a(h)exp|-secd | [a(h)+o(h)] dh 
A h 
iS erin 
h Voces o.+ wh 
te 
h is tke extent cf the atmosphere, which is about akcut 


atm 
30 km, ard by 1s the feight of the laser transmitter. Aisc, 


the wind term has beer expanded to include the relative wind 
velocity due to slewing. Vy» is assumed to be parallel and 
Sevesite in directicr of that ci the target mction. € is the 
angle cf incidence of the wind to the beam so that Voccs(§& 
represents the transverse wind. wh is the effective wind 
generated Ey slewing. wis the angular slew rate. 

Tke £1r0ret pecrtion of equation 2.12 is independent of 
Fath while the second part is not, assuming I(x,y) dces not 
change alcn¢e the prefagation path. This assumftion is valid 


cnly fer a very small amount of L-looming. fhe approach used 


in GUTSAVG is to determine the phase distortion lue to 
thermal blooming by first evaluating the path invariances 
of equation 2.72 . This 1S accomplished fry constricting 
phase screen at the aperture and then removing the hest fit 
tilt, focus curvature, and mean phase. Zernixe polytcnimas. 
are used to model these aberrations. The result is the Tesi- 
dual phase due to thermal blooming alone. The variance of 
the phase is then computed. 

The path dependent term is evaluated within the angle 
interval loop of the program and is artpliceé to tle view. 
ously computed phase variance during each path iteratior. 
The path iteration process 1s diagraumed in the engagement 
geometry section. Also, see Figures 2.1 and 2.2 for a fiow 
diagram of the general treatmert or thermal bloomine ir tne 
program. 

Cnce the total phase variance has been Cetermined, the 
Strehl relation is used to compute the intensity degradatior 


due to thermal blooming. 


= exp(-o7) (2. 12) 
0 


HH 


The resuit of equation 2.13 1S a relative intensity (leew 
ratio. Io is the ideal on-axis irradiance with no phase 
distortion. The Strehl relation above is thought to se umes 
severe a model for fIrel below 0.3 (Ref. 6]. For that 
reason, if o* is less than 1.2 , the Irel will de compare 
uSing polyncmial curve fits developed from GUTSYIF results. 
GUTSMTF is a full wave optics code utilizing fast Fourwves 
transforms. For a description of the curve fit method 
above, see the subroutine BLOOM explanation. [{Ref. 7] 
Combining the thermal bloominy effect with the other 
effects, such aS turbulence] andegitter, 1s accomplished by 


averaging the results of two different approaches. The first 


16 


apprceach is the &SS {roct Sum syuared) method. This rethcd 
Memectiraring tae Prel due to thew™erieccts Of tnermal ticcming 
Mmerne the LIfe gue tc -1tter, Leam juality, and turbulence is 


acccrrliskea as folicas 
i i 
Irel = {1+; ———— -l/| + |—- 1 y 
Tees ( Teer, | peer, (z. 14) 


where irei, Dome@ce ClemllaivcicCGh@ing :esult and Ire]. Loncre 


: c 
result oi tke other ¢fifects. The second approach is a multi- 


pilicative approach are is sim;zly 

Trel = Ir oe 22.15 
The two ccakined Irels cktained Ly these wethods are then 
averaced tc give the total intensity ratic cue to all the 
attenuating or distcrtiny propagation effects. 


tTrel + Lrel 
m Eos 2 ; 


fl a 


= Trel, oy = 


ieee the ideal diffraction Jigited on-axis intensity. the 
Eases fcr the abcve averaginy frecess is empirical in nature 
and is ar attempt to adjust the results obtained by the K&S 
methcc alcne. The resuits produced by RSS were thcught tc 
ke tce optinistic. The multipiicative method, a mcre pessi- 
Mistic aprrcach, was therefcre inciuded. The ultimate [reis 
cbtairedq are very clicse to these obtained by GUITSMIF, the 
full stave ccde. [Ref. 8] 

There are soge linitaticns te the thermal blooming agodel 
used in GUISAVG in addition tc the assumptions already 


menticned. First, the wind (V,) is applied as a corstant 


everyrrtere in the atacsphere. At higher altitudes, tkis is 
rct a gwajer consideraticn due to the higher relative 
velocity cf the Fkeam. At lcw to medium altitudes this cculd 
affect ttexraoal blooming to an extent to warrant the addition 
cf a wind f~refile as a function of altitude. This ceumeaee 
done sitk jJittle effcrt if the data is availaple. Also, the 
Frogram assumes a wind parallel and opposite to the direc- 
tion ci the target metion. This ~reeludes thevcaseuc. 
slewizs¢ *ith the winc and the Creation 2 null see ue. 
Iranscric btlooming, shich viclates one cz the origizal 
assumpricns of the ecvation used, 2S not considered 

[Ref. 9]. For very icw altitude sateliites, the high sler 
rates gererated soulc result in a supersonic relative winré 
across tke ream. In tkis case, the constant [L~ressure a=Ssuiee 
tion is invalid [Ref. 10]. Kinetic cooling ard nolecudan 


Treakccwr are also nct addressed in GUTSAVG. 


Cw. SCATTERING AhD AESORPTICN 


CCISAVE does not contain any type of atmospheric mccel 
with respect to scattering and aksorption data, Extinewnee 
coefficients must be entered by the user in the aprpre;jriate 
subroutines, ALFS and ALFA, or the program can be moditied 
to accért a separate cata file. Once inserted into the 
program, scattering ard abscrption coefficients are treated 
withcut distinction Léetween aerosol and molecular mecha- 
nisms. Therefore, the coefficierts used must represent the 


total effect of scattéring or absorption. 


a(h) = o(h) 9 + off) 7 (2 


1é 


o(h) = o(h) + G6 (ke Cae N85) 


ol aer 


TraneMi¢sion due tc scattering and aosorption are 


computed indentically and given [Ey 


atm 
a 2x (-sec8 sthyan] (2-15) 


Ne 


where 91S the zenitk angle, ard tne inteyration limits are 
meeeiritude on the tlapsmitter and the vertical extent of 

the atnosphere. The extent of the atmosphere in the [Ircgran 
is defined as 30 km. Figure 2.3 Snows toe frojraz aprlica- 


tion ci scattering ard acsorption. 


T. EEFAM CUALITY 


GLISAVG aliows tile user tc describe the pdeam at the 
aperttre in terms of the electromagnetic field anplitude 
feeb ituticr and total ream guality. The intitai agrrlitude 
field used Ly the prcgram is GauSsian in shape with the user 
specifying ¢ waist diameter. The waist diameter is detired 
Ey the lye* point on the distribution. Truncation cf the 
this Gaussian field wall of course depend cn the aperture 
diameter anc the size of the central obscuration. If the 
waist diameter is made large ccu&pared to the aperttre, a 
more vunifora distribution results. 

E¢€am guality at the aperture exit may be specified [y 
two differert parameters. One cf these is the ‘times 
ditfracticn limited rumber', N. N has been used in a 
general way to mean az increase in far-field spot size cr as 
a *power-in-the-kEucket* ratio. In GUTAVG, N is a total Lean 
gGuaiity term. The seccnd beam guality parameter is a nendi- 
mensicnal term representing the RMS phase distortion at tke 
Mecscr wavelength at tke exit aperture, °,.. - 


Soap Gaen 


the parameters are related to €ach ctner ana to the irten- 


Sity dégradation by 


Hy] Hy 


Z 

=. 216 (2.28 

= N2 = 6a Bins 
r 


Ic apply beam guality te the propagaticn rroblen, “ares 
array is ccnstructed representing a phase screen at the 
aperture. ihe phase is assumed to pe a Gaussian randcm vari- 
able witk @ zero mean value. Fer the axi-symmeétric Léarg 
cousicered in GUIAVG, the MTF is [Ref. 11} 


Mi (0) = exe (-Kf*-c, ] } (2.29 
7 1=£ the autocorrelation £unction of Ehe eghace anage 
defined as 
p \? 
Cy (P) = o°exp ae 
) Be | 
2776 
oa BMSe ; 
where o° ig the phase variance, ae 1 and L is “tie 


Fhase ccrrelation length [Ref. 12]. The pean guality MIF 
array is ccmbined with the MIF arrays due to cther prcpaga- 
tion eifects to determine the couplete system MTF and, 
hence, tte total irradiance degradation. Based on the 
Frovided infut parameter, Figure 2.4 shows the general 


treatment cf beam guality within the progran. 


EF. IUCHEULENCE 


GUISaAVG uses the Fufnagel rodel [Ref. 13] for Ce aS var 
indicatcr of the optical effects of turbulence alongeeue 
Fropagation path. Ce is the refractive index structure 


constant and represerts the refractive index in the 


atmcosjthere as a functicn of turbulence induced dénsity 
mlctMaticns. the model is an emfiriagelly derived wertical 
Exotiltle ci oe ; 

Friec [kef. 14] Las develcred a parameter which is 
dicectly related to tke brehavicr cf£ a conerent bead in a 
tuckulent medium. This term is called the effective ccher- 
ence Ciageter, ri: Ir the case of a laser transmitter, [py 
represents a physical limit to the transmitter diameter cia 
near ditfraction-limitéed beam. For a transmitter dianeter, D 
, larcer then ©, , decradaticr cf the nea by turrulerce 
Seemimeceetresrt Db 15 cSiadller than £,05,) nen near ditiracticr- 
limited frcpagation will be achieved. speauw wander or fure 
tilt cccurs for transnitter diameters approximately equal to 
r,- ura (kef. 5] has defined a somewnat different iut 
related tera that car ke thcught of as a lateral cohererce 


lengtt, oy». These tsc guantities are given by 





ame 1.9 , mes (Zee) 
Ci cee peeecs foscn a : 
z 
and 
Po = as esec0 fet dh = (2-24) 
n 
=@— tliat 
a (Zero) 
Oo amg al 
27 
mere kK = Fs | 9 if the zenith angle, and tte linits cf 


integraticn are Dae the vertical extent cof the atumcsfphere, 


and rt ,tte height of the transmitter. 


Po may be specified by the user 4S a program input or 
the program may be ailowed to compute it based 91 ecuation 
2224 . The user may also input 0 ingirectly Dy Specitmama 
the ‘seeing conditions' , a guantity used by astronomers to 
describe angular spread of a steilar point source (See mings 
section and Figure 2.5). Ihe progran uses — to compute the 
atmospheric MIF. 

The NTF of the turbulent atmosphere is deterrired bv 
develoring the structure function of the turbulence. ‘This 
GCevelo, pment is demonstrated by Yura | Ref. 16]. The 


resulting MIF is given by 


_ oma ee. 
M, (0) oS exo |-() | ( ) 


This MTF effectively applies the the ef&rect of ti rbuvemee 
along the propagation path as a phase screen at the 


aperture. 


Fe JITTER 


Beam jitter 1s a user input and 1S specified by the 
variance of the angular excursions of the beam. Jsing the 


7 Sah Na 0 dea ian ieee oie. 
y 22 27) 
20 = 2 eo) ( 
Pp x y 


Winelcy © yeaah oY are the axial variances of the jittered bean 


center in the far-flela= s° sama ey ace random variables with 
Gaussian distribution and in the symmetric case , as consid- 
ered) Dy SGU SA Ce = yy - The resultant intensity distraie 
tion due to jitter will also be a Gaussian distribution with 
on representing the spot radius defined at the 1/e2 
point. In other words, 86.5% of tne beam energy will reside 


within the radius ao .  PpRei ans 


Ze 


Jitter can be shcwn to ke awavefront tilt at the afer- 
Sie misiiewrnse dpe rGCACh, a phase screen at the aperture 
can Le used to characterize the effects of jitter and a 


jitter MIF developed. That MTF is given py 


j==2 2 3S 
M. (0) = ae p (20,,) | (Z. 28) 
8 


Figure 2.€ is a general flow diagram for the treatment cr 


Seeter within the pregran. 


Go. ALAFTIVE OPTICS 


pe doliity €O apEly dd@aptive optics Corrections to tke 
propagation problem has been’included in GUTSAVG. The user 
has several options sith respect to the type and degree of 
compensaticn desired. The fcllcwiny general discussion and 
figure 2.8 frovides the needed insight to the erfects of 
Selecting the adaptive optics cptions. 

Tke wajcr compensation mode provided by the pregrag is 
invoked ky selecting full zcnal adaptive optics with consid- 
€raticr given to isoplanatic effects. When selected, this 
model results in the correction of beam degradation due to 
turbulence. This 1S accomplished by correcting turbulence 
induced tilt and ther adjusting % so as to produce a 
predetergined level cf adaptive optics performance. This 
predetermzined performance is as Measured by the Strehl ratio 
given a residual phase variance determined by the adaptive 
Cptics senscr phase. Parameters determining the phase error 
are tke response pdandwidth of the adaptive optics systeag, 
the number cf system actuaters, the rerlected radiant inten- 
sity of the target, ard the target-to-sensor transmissicn. 
The resultant 9 found in this manner 1s then used tec 
compute the atmospheric MIF. When this ‘corrected! MIF is 
used to compute the far-field intensity, the result will 
represent an adartive cptics ccrrected value. 


Ze 


hithcut invoxiny: full zcnal adaptive Oftics, Siewert 
Hay afpiy a tilt-only correcticn Lor sturbulence tig ae 
case sce cri ali of tke tilt due to turoulence may Le 
removed Lefcre ccmputing the atmospaeric HIF. Tne deyree ct 
tiit ccmpensation is specified by the user. 

As mtéentioned above, isoflarnatic effects are included in 
the acarptive optics calculatiens. This 1S aiso an Ceerene 
however, and isoflanatic calculations may be inhibited ky 
the user. Tke effect considered is the limitation cfr the 
adaptive optics system yivern an isoplanatic angle smaller 
than tke target lJead angle. ~hrved [ket] ie a eoy a 
discussicn cf iscpianatism and develocment cr the iscpia= 
Datic eabcue: 

Altkcugi adafctive optics cecmpensaticn Lon eis 
kiocrinsg is not modeled in a strict sense, a therzral 
riocBing ccrtrection factor can te applied. Inis ftacter mie 
Simply a fractional ccnstant that multiplies the thermam 
Elooming phase variance befcre the Strehl relation is used 
to ccupute the intensity degracation. 

Refer tc the input definiticns and the sSucroutines 
involved with adaptive optics. In addition, refer tc figtre 


ze/ £cr wore explanaticno. 
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Fagure 2.9 GUISAVG Thermal Blooming Algorithz. 
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Figure eae GUISAVG Scattering and Absorption Algorithn. 
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Figure 2.4 CUISAVG EFeam Quality Algorithg. 
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Figure 2.5 GUTSAVG Turbulence Algorithn. 
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INPUT DECISIONS 
NELAGA FULL ZONAL AO 


1 YES 
1) NO 
NOTSO INHIBIT ISOPLANATIC CALCULATIONS 
l YES 
0 NO 


ADAP RESIDUAL TURBULENCE TILT 
0.0 REMOVE 100% OF THE TURBULENCE TILT 
1.0 REMOVE OS OF THE TURBULENCE TILT 
AOBLOM THERMAL BLOOMING CORRECTION 
0.0 REMOVE 100% OF THE THERMAL 8&1 OOMING 
1.0 REMOVE O% OF THE THERMAL BLOOMING 
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Parque?) Za 7 GUISAVG Adaptive optics Algorithn. 
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Figure 238 GUISAVG Adaptive Optics Algorithm (ccnt). 
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GUTSAVG Adaptive Optics Algorithm (ccnt). 


Figure zond 


TLTISO 


COMPUTE WAVEFRONT TILT 
DUE TO ISOPLANATIC 
EFFECT 


YES NFLGSN= |} 
NOISO=0 





SINTLS 


COMPUTE IRRADIANCE LOSS 
DUE TO SCINTILLATION 
EFFECTS FOR AO 
APPLICATION 














YES NFLAGA= 


ML GSN#-0 


< 


MTF OIT 


RECOMPUTE JITTER MTF TO 
INCLUOE ISOPLANATIC TILT 


NO 








GUTSAVG Adartive Optics Algorithm (ccnt). 
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Figure z.11 GUTSAVG Adaftive Optics Algoritha (ccnt). 
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Figure Zane GUTSAVG Adarftive Optics Algorithm (ccnt). 
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III. PROGRAM AND SUBEROGRAH DESCRIPTION 


Pee f£BCGEAS INPUIS 


The following are the inputs necessary to utilize 
GUTSAVG. Default valtes are iandicated for parameters not 
absclutely rejyuired fcr program operation. At NPS these 
input farameters are entered via an input rile. A cCcEYy oF 
mags file if provided in Appendix As {*) indicates a ncndi- 


wensicral f;arameter. 


2 [IA {meters) 


Thre Ciaazeter of the trarsmitter aperture. 


e [TIACES (meters) 


Tte Cianyeter of the central obscuration of the trans- 


Hitter. 


« EFAMSZ (meters) 


The Gaussian waist diameter of the amplitude distritu- 


ticn at the aperttre. (measured at the 1/e2? point) 


e WAVE (meters) 


Tre laser wavelercth. 


a) FICIAL (Watts) 


Tre tctal power at the aperture. 


e TLFIMTI (¥*) 


This is the c£ter used “tities tdat teraction pec 
nurker® and represents total ream guality. As used in 
GCIS, it is related to the RMS phase distcrtion at the 


aperture by the Strehl approximation 


i 276 ° 
(TDFLMT) a exo |-( a | (3.1) 


(IDFIMNI)*2 , therefore, is equivalent to the tcatic Cima 
crmaxis diffraction limited intensity to the on=axie 
intersity resulting from the near-rield fnoase distcr- 
SGU a oer 


ee: 3 
(TREE) — en, (sag) 


Fer @€ discussion cr the ligitation cf 3.1, see 
[bef. 1S}. 


« WAVEEC (*) 


This term is the EMS phase distortion at the laser aper- 


ttre, nendimensicralized Ly the wavelength. 


5 se 
WAVEBO = = oe 
e SCALEC (meters) default = DIA/S5 


This térm is the transverse phase correlation lenctt at 


tke aperture. 


ae 


Serio be (acCcsec) default = £(RHOOQ) 


This is a qualitative term used by astronome:zs to 
describe ‘seeing conditions* in the visible range. If a 
point source is viewed from the earth, it may not appear 
aS a point source but as a 'smear' or spot. The angular 
spread Of this Spot is the parameter THSEE. Tf RHOO is 
Mewespecitied aS an input and THSEE is, the program will 
lisemitorheto sSOompute RHOO. Other than this case, THSEE 
is not used but will be computed as a function of RHOO 


and included as an output. 


e HGRND (meters) 


The height of the ground at tne transmitter position 


above MSL. 


@ HTRANS (meters) 


The height of the transmitter above. 


e HSAT (meters) 


The orbital altitude of the satellite above 4SL (at 


Zenith) « 


e THETMX (degrees) 


THETMX is the angle measured from zenith below which the 
laser will not transmit for a zenith pass. For offset 
flight paths, it should be noted that the transmitter 
will point below this value. For an illustration of the 


engagement envelofe, see figure 3.1 


oro 


e iC iE (meters) 


Fer @ target that dceS ret fass direct overhead, this 

irzjut specifies tke amount that the target gqrecuncweeaes 
is cfiset frem tre ground track of the overhead case. It 
is tke distance, as measured from tne transmitter, toa 


Ferpendiculaz intelsecticn cn tne Ggeoutd trac. 


@ KHCO (meters) default = £/ SP!) on 2507. 


inis is the turbulence coherence length as defined Ly H. 
- Yura. [Ref. 720}. (see turbuLence section fcr ere 
oS 


a : 
discussicn ard references) 


e VO {n/SEc) 


Tke atacspheric sind. The direction or the wind is 
parallel and oppcesite tc the direction of target ctageee 
Ncte, in the present yprcyran, VO is a COlsteant 2a0ee 
G€nt Cf dltacude: 


* ACELCM (*) 0.0 to 1.0 


This crarameter allews the user to correct for thermal 
Llecnansg as if by adaptive cprtics. The value Entered say 
range from 0.0 te 1.0. If 0.0 is used,comprlete therzal 
blecairng compensation will cccur. Conversely, if 1.0 is 
ertered, no compensation will te applied. The variance 
cf tke f~hase distcrtion is tuitiplied by this ccrreéecticn 
factcr Lerore the Strehl] relation is used to compute the 


relative intensity reducticn due to thermal blooming. 


= = exp (-0* * AOBLOM) (3-4) 
0 
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¢ AVGSET (*) 0.0 to 1.0 


AVGSEIT allows the user latitude in defining tne tar- 
field sfot size te other than that indicated by GUTSAVG 
analysis. A value of 1.0 would result in the peak irra- 
diance according to the progran analysis. An entered 
Varwe Of «9, LOE example, would resuit in a peak irradi- 
ance 504 less than the program would otherwise indicate. 
AVGSET, then, iS an adjustment factor that allows the 
user tc account fer erfects not addressed in the GUUZSAVS 


prepagation calculations. 


«e SIGJII (radians) 


PeGull is the on Var tancourOnesommrerug and tracking 


wie ter. 


¢ ADAP (*) 0.0 to 1.0 


mimes telm fs aGccEEeCCtl Gn Lactor for tilt due toe turbu- 
lence. It represents the residual tilt after AO conpen- 
Saticn. ADAP may be varied from 0.0 to 1.0. A value cf 
1.0 would result in no tilt due to turbuience pDeing 
remcved while a value of 0.0 would result in total tilt 
due to turbulence ‘compensation. If full AO is selected 


the ~rcegram will set ADAP equal to 0.0. 


e NFLAGA (*) O cr 1 


NFLAGA is a Selection indicator for full zonal adaptive 
Crtics. Enter a Q if AC is not desired, enter a 1 if AO 
1s désired. If AC is used, XJT,BWIDTH, and NA must ce 
specified. 
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2 NCISC (*) OVcr 3 


ACISC is a selection indicator Zor isoplanatic cate. 
ticns. Enter a 1 to anhibit asoplanatic calcu lage 


enter a 0 if iscflanatic calculations are desired. 


« XJ1 (Watts/<sterad) 


XJIT is the target radiact intensity. it is one of the 
factcrs used to céetermine how noise arifects the restore 


cf tre 40 system. 


* EWIDTE (Hz) 


ERILTIH is the bancwidta of the adaptive optics systen. 


e NA (*) 


NA as the number cf AO actyuaters used tc ferticrm phase 
a 


cjvstnrent. 


*« AESICZ (%*) 


This is the percent transmission at zenith at the 
sensing wavelengtk of the AO system. This term is used 


in tke deternination of retadtive noise at the AO Ser. 


° NT (*) 


The rumker of iteration steps for the angle lecr cf tte 


Prcgraai irom THEIYX too: 


* NZ {*) 


Peo mvUlsecheOcsa titude IDtervals (fOr alsorpti or and 


SGQGterIngydeteriiratior. 


* N3 (*) 


ree runier, Of altitude antervals £or tae turbulence 


@alcutations. 


~ NY (*) 


TkKe number of iteration intervals ior tne ATF calcula- 


Eaecnse. 


e NS (*) 


Tne nugker of intervais fcr the slant path update cf 


thergsal blooming. 


Ej. ERGACEMENT GEOMETBY 


The target engagenent gecmetry is that of a earth-Lased 
transmitter and a taroet Satellite at a given orbit. Ne 
attempt is nade to deiine an engagement envelope Lased cn 
thresrold irradiance cr fluence!. The input parameters 
defining the engagemert window in GUITSAVG are THETMX, LCF, 
and HSAT. Tre general gecmetry of this window is shcwn in 
Figure 3.1 . Only half of the total transit window is 


addressed in the progragz caiculations; the results are the 


i The GUTSFP (fcctprint) version of gtts was writter 
Pomac tiisam™ 3XCepe fCr this feature, the Propagation calcu- 
jdaticrs are identicai to GUISAVG. 


4 


dad 


same fcr €ither nalf. The prograi owtput f£ereect oem 
halI~window evaivaticr except for parameters such as tctal 
fluence and Shot time which are simply double the cciguted 
values. 

Figure 2.2 defines the earth center angle (ECANG) nea 
is a functicn of the tser input THETMX. Most of the gecre= 
tric calculaticns are referenced to earth center. Therefore, 
this argle is used £cr computin, such positioning data as 
the angle irterval at which the irradiance will ve evaluated 
(see Figure 325). 

Cifset [Cliyat paths require a coordinate transiatica as 
Snown in Fisure 3.4 . Positicn ang velocity relative Coma. 
transmitter are comptted as in Figure 3.5 . It skould Ee 
noted that if the flicht path is offset, the zenith angle 
Will exceed THETMX fcr fart or all of the window. Inveeee 


tecause tre windcw 1s Gefined in the x-z plane only. 


qu 


He - Zenith altitude of target 


= - Maximum zenith angle for engagement (Zenith path) 


L. - Offset distance 


Offset 


é a 
Zenith E likeflae Pehelel 


mibroght Path 


Pegagement 
Window 


Se 


Figure 3.1 General Engagement Geometry. 
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x 
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= cos 2 sing. 9 5 
S 
R Se 
S - Target satellite 
T - Laser transmitter 
a - Earth center angle 
us - Maximum zenith angle for 
engagenecn 
Ro - Earth center-to-target distance 
R - Barth center-to-transmitter 
distance 
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Paguresss Earth Center Angle (ECANG). 
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Itracdiance is Computed at 
interval midpoints 





The number of intervals 1s 
a iscrmenpl eye moe Le 


Pagure. 3.3 Angle Intervals for Irradiance Evaluaticr. 
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Offset is from T ems 
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Figure 3.4 Coordinate Translation due to Flight Path Cffset. 
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Velocity components 


Va) eos (6. 
x Ss e 
; ae = Vsin(6_)cos(y) 
3 = Vsin(6.)sin(y) 
y 
x 
S Position 
Ke= Ra sim(s.} 
S e 
Y = Ly-cos(y)+(4Z)-R-Lz) 
-sin(y) 
Z = -Ly:sin(y)+(Z )-R-Lz) 
"GOs (y ) 


where 2Z,= R_cos(é@_) 
Ss e 


Range from transmitter(T) to 
Sacre lite (Ss) 


= ag a2 
Ri Z Xx % 





(For Ly and Lz definition, see Figure 3.4) 


Figure 3.5 Velccity and Position Relative to Transmitter. 
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C. MAIN FECGRAM FLOW DIAGRAM 


Tké€ waiL proyram gestly consists of geometry calcula- 
tions and decision ficw poircts. The decision points allew 
ktranchirs te adaptive optics and isoplanatic subroutine 
Most cf the propagaticn calculations are done withizn yop 
Ssubrcettaines. Only tae major iteration loops reside in the 
Bain freocran, the angle interval loop and the combined WiF 
integraticr loop. Tte foilcwinrg is a yenetal £104 diag 
for tke naiz program. Decision variabies are user inrputs cr 


Frograg cefined fparaueters. 
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ETgune 3. 6 





RE QUIRED 
PROGRAM 
INPUTS 







COMPUTE EARTH CENTER 
OFFSET ANGLE ANO 
COORDINATE TRANSLATION 
DUE 10 FLIGHT PATH 
OF FSET. 







COMPUTE ORBITAL 
VELOCITY OF TARGET 


COMPUTE EARTH CENTER 
ANGULAR RATE OF TARGET 


COMPUTE TARGET 
VELOCITY RELATIVE TO 
TRANSMITTER 


BASED ON THE TMX ANO 
ANGULAR VELOCITY 
COMPUTE SHOT TIME AND 
EARTH CENTER 


ANGLE LIMITS 


UOCST ; 


DETERMINE THE CONSTANT 
THAT CAUSES THE 
INTENSITY INTEGRATION 
TO RESULT IN 


THE SPECIFIED PYOTAL 






Main Prcgram Flow Diagran. 


Se 


FARFLD 





COMPUTE THE 
FAR-FIELD IRRAOIANCE 
AND THE 
CORRE SPONO ING 
APERTURE MIF 








WAVE BQ=0 





DETERMINE THE WAVEBQ 
THAT WILL RESULT IN 
THE SPECIFIED TOFLMT 


BOIREL 


COMPUTE TDOFLMT USING 
THE SPECIFIED WAVEBQ 



















CONSTRUCT THE BEAM 
QUALITY MTF ARRAY AT 
THE APERTURE 






CONSTRUCT THE JITTER 
MTF ARRAY BASED ON 
THE SPECIFIED JITTER 








Figure 3.7 Main Program Flow Diagram (coat). 
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Ww) 


ABSORB/SCAT 


COMPUTE THE 
ABSORPTION AND 
SCATTERING FOR 


ZENITH TRANSMISSION 





RHOO0<0 





RHOTRB 


COMPUTE THE 
TRANSVERSE COHERENCE 
LENGTH 





NFL AGA21 
NOTSO=0 





YES 


ISOTRB 


CALCULATE THE ZENITH 
ISOPLANATIC ANGLE 






Figure 3.8 Main Pregram Flow Diagram (cont) ~- 
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Figure 3.9 


COMPUTE THE LOG 
AMP IL ITUOE 
SCINTILLATION FOR 
ZENITH TRANSMISSION 


THIS IS USEO LATER FOR 
AO CALCULATIONS 








PRTPHS 


COMPUTE THE THERMAL 
BLOOMING PHASE CHANGE 
INDEPENOENT OF PATH 






RE SPHS 


COMPUTE RESIDUAL PHASE 
DUE TO THERMAL 
BLOOMING BY REMOVING 
MEAN, TILT, ANO 


FOCUS CURVATURES 











COMPUTE THERMAL 
BLOOMING PHASE VARIANCE 





BEGIN ANGLE ITERATION 
LOOP FOR TARGET PATH 


INITIALIZE ACCUMULATED 
FLUENCE 





FLUX=0 


£4 


PROVIDE 
ZERNIKECOEFFICIENTS 
FOR POLYNOMIAL FIT 


OF MEAN, TILT, ANO 
FOCUS 





Main Program Flow Diagram (cont). 


Figure 3.10 


COMPUTE ELAPSED SHOT 
TIME AND COMPUTE NEW 
EARTH CENTER SHOOT 
ANGLE 


COMPUTE THE TARGET 


COORDINATES 
REFERENCED TO THE 
TRANSMITTER 


COMPUTE RANGE TO 
TARGET ANO 
INSTANTANEOUS SLEW 
RATE 





COMPUTE ZENITH ANGLE 
AND ANGLE OF 
INCIDENCE OF WIND ON 
BEAM 







AV 





COMPLETE THERMAL 
OOM ING 


BL 
CALCULATIONS 8Y 
INCLUOING PATH 
DEPENDENT EFFECTS 







Hain Pregram Flow Diagram (cont). 
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APPLY THERMAL BLOOMING 
AO CORRECTION FACTOR 
AND COMPUTE THE Ire} 

DUE TO BLOOMING 








IF THERMAL BLOOMING 
PHASE VARIANCE IS 
>1.2,CURVE FIT TO 
GUTSMTF RESULTS TO 


OETERMINE Irel OLE TO 
BLOOMING 










ADJUST SCATTERING, 
ABSORBTION AND RHOO 
FOR SLANT PATH 






COMPUTE THE LOOK-AHEAO 
ANGLE 





Figure 3.11 Fain Pregram Flow Diagram (cont). 


NO 


YES 


TRBREL 


DETERMINE THE Irel FOR 
A SYSTEM CORRECTED FOR 
TILT ONLY 









Ire? = EXOTST 


COMPUTE THE Irel OF A 
PERFECT AO SYSTEM 
SUBJECT TO NOISE 


Irel = EXD 





YES 


SET THE NOISE FLAG 
TO INDICATED SYSTEM 


ONL Y 





Figure 3.12 ' Nain Program Flow Diagram (cont). 
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COMPUTE ISOPLANATIC 

ANGLE ANO TEST, BASED 

ON LOOK-AHEAOD ANGLE, NISO=1 YES 
TO DETERMINE IF 


FULL ISOPLANATIC NISO*0 NO 
CALCULATIONS NEED TO 
BE DONE 





YES 










NISO = 0 





COMPUTE IRRADIANCE 
REOUCTION OUE TO SLANT 
RANGE. 


YES NE LGSN= 


NO1SO=0 


VES 


TLTISO 


COMPUTE WAVEFRONT TILT 
OVE TO JSOPLANATIC 
EFFECT 


Figure 3.716 Main Program Flow Diagram (cont). 





SINTLS 


COMPUTE IRRADIANCE LOSS 
DUE TO SCINTILLATION 
EFFECTS FOR AO 
APPLICATION 






MTF ATM 


COMPUTE Tec eee’ 








MTF JIT 


RECOMPUTE JITTER MTF TO 
INCLUDE ISOPLANATIC TILT 





Figure 3.14 Main Program Flow Diagram (cont). 









COMPUTE THE ISOPLANATIC 
MTF FOR THE FULL AO 
APPLICATION 


INTEGRATE JITTER, BEAM 
QUALITY. TURBULENCE, 
APERTURE ANO 
ISOPLANATIC MIFS TO 


FINO Irel 










CONTINUE THE MIF 
INTEGRATION LOOP FOR 
0.0 < xX < DIA 


INTEGRATE JITTER, BEAM 
QUALITY. TURBULENCE 
AND APERTURE MTFs TO 

FINO Irel 


Figure 3.15 Main Program Flow Diagram (cont). 
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COMBINE BLOOMING Irel 
ANO Irel FROM OTHER 
EFFECTS BY RSS ANDO 

MULTIPLICATIVE 
METHOD 





RELTOT 


AVERAGE THE TWO TOTAL 
Irel TERMS FOUND ABOVE 


CONTINUE ZENITH ANGLE 
ITERATION FOR O<THETA<THE TMX 
APPLY SPOT AVERAGING 
FACTOR IF SPECIFIED 


THE NUMBER OF INTERVALS = Nl 


DECEASE IRRADIANCE 
BY SCINTILLATION 
LOSS IF FULL AO USED 


ACCUMULATE FLUENCE 


COMPUTE ZENITH ANGLE 


COMPUTE TOTAL 
ILLUMINATION TIME 


CONSTRUCT QUPUT ARRAYS 





Figure 3.16 Kain Program Flow Diagram (cont). 
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Eee CEnCUTINE DESCRIETICN 


femeaecCns and SCAT 


As noted in tké preceding chapter, absorption and 
Seatterinc efrects aze treated indentically. The AESCEE ard 


SCAT algcrithas ferfcim the integration 


h 
atm 3.5 
T, = exp ai a(h) ah i>) 
De 
eee 
T = exp a Gena (3.9) 
S hy 


Tne rcesult is the trarsmission due to the totai akser;ytion 
Smmtctal Sscattezing. The correcticn foe slant ranye is 
@pplied sithin the angle interval looo of the main frograiz 
and is siafply 


(sec8) 
of = 7 (3.7) 


is the zérith angle of the target with respect to the 
teansaitter. 
fhe extiancticn coefficients for absorption (a) anc 
for scattering (9) are provided to the subrcutines Ly a call 
tO Ecttines ALFA and ALrS respectively. 
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mms See em a i 


| 
D Bee eee - 
| AESCRB and SCAT Progragz Variable Definitions 


intervals 


( 
| 
i 
Vériakle Fortran nage eran: | 
Sub Warn Definiticn 
l I_ (AESCEB) T TABSC Total golecular | 
| A transmission 
i Ts (SCAT) z TscA re Total scattered { 
{ transmission { 
{ am | 
| h {feters) HATMC HATIMC Heignt of j 
j atm atmosphere 
| i, (meters) HT HTRANS deight akove MSI 
{ OF transmitter 
c (b) eons AL? * ADSOELtICR (Coem | 
| at given altiaecue. | 
i et | P 
j go (Ek) (ki) ALS: = scattering coe: { 
| at given altitude | 
i . . 
= N N2 # OL Integ¢ratcy-. | 
i { 
4 


nn 


These twe rortines are also identical. They frevide 
the extirction coefficients tc ABSORB and SCAT for a Speer 
fied altitude. At eack altitude where a coefficient valte is 
desired, a linear intéerpolaticn is performed between data 
FOlntS sup, lied fy ere user. Data statements precede each of 
these sukprcgrams, ard it i¢ with these statements that the 
abscryticn and scattering data should be entered. Urits fcr 
the ccefficient and ccrresponding altitude should be km and 
km respectively. 
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| 


TADCE eis 
ALFS and ALFA Eregram Variables and Definitions 


cr peel 
| | 
i 
! ! 
| ee 
{ Variakie Fortran hae ee | 
DUE Nain Definition 
O(k) {key Ss Ages SGaditeuc lanigmsc oe £ < 
| at specifiea | 
altitude 
j o (bE) (ka): A ALE MDSOretio0n ¢coce . | 
| ceeoweeit Lod 
aititude 
{ h (metecs) H = Seecwr1c alritude { 
| 
— ALT(XL) - Altitude data | 
j dist 
_ RSA 5) AOsOmeelLon data j 
List | 
j Sar AT Sah) = Scaiteking data j 
{ Bins fe | 
i 
| a NL = Number of points { 
Leena ats te 

| 
| | 
Be cee SRE ee a a a a | 

-- £0Ccsi 


Fer a gausSiar béam with constant phase, the ititial 


amplitude distriktuticn at the transmitter aperture is 
WG) = voexp |-(e/w) (3-&) 


where U.is the amplittde and wis the spot size ] Ref. 21]. 
Maewceurecse of J0OCST as to ccmpute the Constant U, for a 
given aperture power. 

Ey the scalar wave approximation, the intersity 
Srotriktuticr 1s [Ref. 22] 
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Py eee (3.9) 


fo relate the field” distributicn to whe power, the intetogme 


distrikuticr 1S intecrated over toe apertiuzce. 


‘s 
. 2 
Py = 6 fra: 3.10) 


Ee 
i 


Substituting 3.8 @nd 2.9 dan tom ead 


ae 
© 
are 77 U4 exp | ~ (2/1 ‘ ita (2547) 


cr. 
a 


then Leéalranging teris, produces tne expressicn for U.. 


ane aes 

= Py rf xe Lex} aire (3. tee 
ioe 
i 


The iantecration limits Eo and tI. are the radius of tke 
transnitter and the ckscuraticr, resnectively. 

COCST perforus the integration in 3.12 using the 
trapezoical rule. A call is made to Suoroutine FIELLE tc 
evaluate exp([-(ryvw)?] which 1s simply U(r) with U,= 1. For 
this reascn, U,is defined as unity when FIELD is called [Ly 
VOESi. 
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mm ma mm mr ce ee se mmc mmm mmm ree cm em i 


TAB read 


DOCST Frogram Variakies and Definitions 


om oe @ ome Pin pen 





; 
| 
j 
i 
{ 
j tee pe 
{ Vsriari¢ Fortian fale - 
Sue "Hain Definitica 
| beh , U0 = Normalization | 
j GON scant | 
| j 
| Fe (hkatts) P PTCTAL Total aperture i 
| exit power i 
{ —" | 
| — N ret moe 1 Cera tao 1 | 
{ increments 
: i 
r_{reters) nO ~ Orem racLus ice 
j O transi tter j 
_ i i 
| rj, (meters) KE = aden US 6 T j 
} o opscuration i 
j 7 = : - 
U(r) (nw) a Uk = rie€elc ampiit ude | 
at cadius r j 
d 
w (meters) BARAL - Kaas Of Sp ot } 
| dtwaper ture | 
{ 
| 
| | 
L a a es es ‘cence ees ams ae ces ees Sees ees Geshu soem cemat-aii ais suai seunomsemans pine ad 
4. FIELD 


FIEFID caiculates the field distributicn for a axi- 


symmetric Gaussian beam wits constant phase. 


U(r) = oexp |-(r/s)* 3-13) 
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Ca: See er ee eee 


TABLE IV 


FIELD Prcgram Variable Definitions 





| 
| 


i 
t 
i 
| 
j 
| : 
1 Nariatle Fortran hame _—. 
; Sub Main Definition 
i U(r) (i) mn? UR = Fiea anplieeaee 
| at Cc i 
bs 
U, (hk) me UO - Normalization 
| constant | 
| w (meters) BMEREL - Radius of spot { 
| . at aperture j 
| j 
{ pn 
Sc. FA 


In the Fraunhcfer réegicn, the agplitude distrikutici 


can Le icund by takinc the Fourier transform cf tne apertrre 


cistrikuticc [Ref. 235}. Using pelar coordinates and nota 


axial symmetry, U(p) in the far-field can be expressed as 


= = f foe exp |- (#32). O- cose] rarae (2. Tas 


tsing integral relaticn for the zero order Bessel function, 


- 


egquaticn 3.14 can be simplified to [Ref. 24] 


U(o) = 2 fy U (as) 54 (2324) fae {2. Ty 
a 
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Piemerrtclcanyeumetme Llarsticld 1s then given ty 1 {r)={U{r)!¢ 


Gein tErns of ecuaticn 3.15 


> S 2 
I(o) = (=) S ote) a (2228) cae | Em 16) 
a 


mmere £ ang a are the radius ci the trarspitter and the 





Tadius cf the obscuration, resfectively. [Ref. 25] 

FARFLD evaluates 3.16 at a specified number cf 
increments in the far-field and assigns these values tc an 
array F(a). F(a) 1S [crmalized with the total aperture fower 
Pomerat the MIF ~rodvtced would be unity at the origir. Acte, 
however, thet when tle MTF is ccmputed in this case, the MIF 
is nermalized so as tc frodtce an Irei (intensity relative) 


value wren inteyratec. 


ie 
E 


ee 2 te — ae re aa S 
M_{?) ers ve 1D) Io (27S odo (2.17) 


ee 


Sees Ceternaned by evaluating equation 3.176 with r=0.0. 
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€. ECIBEL 


TABLE V 


| 

| Weed are Fortran name 

| Sub - Main 
em CKR - 

| 

l p R1 - 

| 

j E RO a 

| 1, {katts) PMAXCG PMAXO 
j 

| I(x) (Wym2) F(T) TEESE 

| M (Dp) G (I) IREMIF 
| - Dx Dx 

| 

l F,(hatts) P PTOTAL 
i 


FARFLD Program Variakles and Definitioms 


Radi Usain 
faceeaeld 

Radius ia 
aperture 
Far-field on-axis 
pies speck re te (Oye al 


di tiractziom 
limited Lean 


Irradiance in 
falazielid 


Aperture MIF 


ApenCcure MiG 
increments 


Total power i1n 
the aperture 


See Min OE tees Sie tte EEK Gees Meee ete Ge tte, De eee OE Ge Bo a, i ee ee on Bee fe, Oo eR ee 


| 


If the "times diffraction limited number' (IDFIMI) is 


not specified, this reutine computes the Irel value due tc 


Feam guality which is then used to compute TDFLMT in the 


Bain f-regram. The calculaticn performed is 


igs = 2 fue My 


where Has the aperttre sian aud is the MTF for the ckean 


quality ~hace Setecn sie 


pdp 


(3. 18) 


is defined in subroutine FARFID so 


that this integration froduces an Irel value- HM) is thewecua 
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g 


Poets IF aS ccmputed es, sukrcutine HIFBQ. IDFLMT is 


Semeuted 1n the Nalin ~rogram and is 1/ Irel . 


pe ee i ls dn i _eagn  aeee ei Ga Rei ele 


le ee a fee oo — . .—.  .- ai 
| 
TABLE Vi | 
ECIREL Precram Variakles and Definitioms 
Vartiarle rortran nhawe ne 
Sus Waza Definition | 
N = Ei ee Times diff. { 
iigited # { 
| 
Irel REL TEC = 
¥) er) ILRRMIF Aperture MT? | 
= | . 
My (P) Be) as Whee OL Ehase | 
screen 
= N4 N4 # Of iteraticns ~ 
POmemer CAalculaticns 
{ 
i 
ee ee ne ee nn ee a SN Sa SN Ma Se SE SS LS A ED A <a we wD es we eee nl 
7. LEDKAY 
If the RMS phase distortion parameter (WAVEEC) is 
not specified, this reutine ccmputes it based on the 
apprczimaticn 
Trel = exp (-o07) BAD) 
276 ns 
where oO = r >and 6 emg 1S the RAS value ci the phase 
distcrticn at the aperture. Irel 1s the intensity degrada- 


t 
oT 


c 


ion due to beam guality and is egual to 1/(N)2. N is the 
mout farameter IDFLM1, DETIWAV first evaluates o° using 


guaticn 3.19 and ther uses this value as a starting feint 
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to ccapute a more acctrate o* by iteration using sukrcvtine 


Je peek o* as used tc Compute the beam Jjuality MIF and tke 
keam cuality Irel aS in eyguaticn 3.15 . 6% is adjusted yuan 
the frocess repeated untii the Irel found by equaticn 3.78 


is egual tc that determined Ey 1/(N)?2. 


ae a eS a ee 8 RS Se Se a ee Se Se sa eS ee SS ee —" ae 


TAELE Vil 


DETWAV Program YVariatles and Definitions 





Veriatle Fortran wane oe 
Ine 2 u ~ Hain Perinition 
——— VAREC VAREC Constant | 
6 WAVEEQ aAAVEEQ Phase | 
_ims distortion | 
r parameter j 
1/N2 Ree 0 - 1/ (IDFLMT) 2 ‘ 
| 
{ 
ee 


P= (iii term mm Mc, MMi i Ty ih cpanel om a Hn Es a! Ms ee 
[ 
| 
[ 
| 
| 
\ 
| 
| 
| 
L 
i 
| 


E. FHIFEQ and Bm 


These two routines calculate the beam guality MIF 


array. Tie £2F as f fet. §267) 
M,, (P) = exp (-K* [o" = c, ))) { a2 


where ale is the autocorrelation function of the fhkase and 
is tke phase variance. C (0) igs assumed to-be Gaussian. 


lettirg Cy(o) = o*exp[-(0/L) 2] , where L is the phase 


correlation length, results in [Ref. 27] 


val 


270 ms Sa 24 
M, (BD) = exp ee 1-exp[-(/L) J (J-<1) 
° ems 
i) PemtiewwayClensinenmis phase GiStortion and is ar 
Gpput farameter (WAVEEQG). IL ,also a user input (SCALEC), 
will default to 1/5 tke diameter of the aperture if cct 


ctherkise specified. 


seams rae 
4 
| TAELE Vili { 
. ee i 

MIFEC and BM Erogram Yariables and Definitiors 
i 4 
i Variakis roOrtran fare . 
i, 2 aus Hain Definition 
T d 
| rms VAREC VAREC Bean juality 
K variance | 
¥ (9 Eon) BC eee Beam gualit: | 
| pl ° ceed u ! 
| - N4 NY ao 
| Pree UraciOR | 
ancecrenents 

| k (m) ° = - (27/2 } 

| 
j { 
4 a ee ee Sp lace Same ee RRO il a et See aoe me ames 





v2 


ISCTIRB calculates tne zenitn isdplianatic argle Gem 
use in the adaptive cjtics portion of the progran. The ancle 


iS Given Ey [Ret s2ean 


atm ; 
- 2m ; Stee | 
Oy =e [$25 ccs fog h 2 


as 


(2.28 


For the rear zenitn case, secod=l1. Ce is the refractive 
gndex structure constant and is calculated ry subroutine 
CN2H as function of altitude. See fricd)] Ret. 295 eico 


discussicn cf the abcve angle. 


Ee Fg 
| i 
TABLE ix 
| ISOTRB Precram Variatles and Definitioss | 
i er | 
1 WVariarles Fottran Name Pe 
| Sub Wain Definition 
Og (rad) ISOANKE ISCANO 2enith 
iSO plaice 
j amc Le j 
\ ~ 2, ; { 
C2 (ng) ° CN2 _ Refractive index 
structure 
l constant { 
| . l 
| - N N3 Inteyration > | 
{ intervals foc | 
turbulence 
| . { 
@ (dec) = OMEGA Z¢nith angle | 
| | 
a i es ee a a Ss <a eae ie =O Fu ee 


10. SCINT and SiNILS 


if full adaptive optics are used, amplitude scintil= 


laticn e¢eifects will act to degrade the perfcrm@ance of tke 


ve 


4 
Ce 
}- 
2) 


conventicnal phase ccrrectirg adaptive optics system. 


Subrcutire computes tke log aaplitude variance or the 


Cc in 
in a 
@  b- 
oO mf 


Biblatices LOL Zenith transtissicn. The approximaticn 


is 


Ge = 56k 
Zz 





For tke cfif zenith case 


o? sec(8) 
Z 


N 


wnere 90 is the Zenitk angle {Rei. 30]. Ine relative ifraci- 


mmee 1CSs Gte te SCirtiliation is 


oe, 
as 
a 
A) 
cs} 
ae? 


AMP LOSS = exp (-o3sec() 


Miemccrrection £ecr tle orf zenith case is applied it sutrcu- 
tine SINILS in the ancle loop cf the main program when full 
adaptive ocptics are specified. This loss factor has keen 
Shown tc te limited tc approximately .55, thereiore, a 
Gefault value is used in the air program in tne case where 
Sepmeeacnr 2.25 yields a value less than 0.5 [Ref. 31]. 
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SCINILS and SCINT Pr 


fee ee 


We tele tele Fortrab 
aus 
kK (DD) CK 
| A 
| C2 (fm) CN2E 
| 
l 
| a 
] fe} siG %Z 
] x 
Agr lcss TAME 
j 
{ 
tL 
11. CNZE 


mmm ne se re ee re ee ee em ee ee ee 


TABLE 


cgram variables and Defimiticns 


ale Pe ae 

Haiv Deti nit ren 

= 2 a 

=< Refractive index 


Structure 
constant 


SIGX2Z Log amplituie 
Va iaeion 


TAMF LOES Caused [Ly 
SCintivilatrvon 


CN2EH computes the igndgem of FecractiOnest: vere 


constant as a Euncticn O£ a lei tide. C2 is determined using 


Hufnacei's model witk an added term to include the near 


surface effects. [Ref. 


Ore me | Sal) = 
Nn . 


oo 


32] [Ref. 33] 


-h 
yh exp (585) (3. 26) 


(27a )exp (50 


+ (10 esp (58 


Here, the urits of h are meters. 


1z. EHOTRE 


This subroutine 


computes Yura's [Ref. 34] lateral 


coherence length,Po , in terms of Fried'*s [Ref. 35] conere 


ence Cianeter, Co. 


fee 


ee a a SU Oe 





Ce OR es Be Be gg, ee Ge ee Be ne ee Oe poe! 


Ae eeee Ree gae ea 


Dg = eo (se) 


Em@gstatutincg Fried‘s definiticn for cr 
Q 


—3, 
1 Se en : (5 
= tea Z 
ain. Gmc. E C,, (A) dh 
ts 


tar 


Nd 
OD 
we, 


| 


Baer aa 


RHOTRB Procram Variatles and Definitiozs 


Varlerie Forirau sane ae 
Sul Main Definitics 
Keeti)s + CK = 27/7 A 
—2, 

C2 (nr) ° CN2 - Rerractive index 

n structure 
constant 

O, (rete cs) RHC RECO Conerence itengtk 

= N N3 t (Om 1LNtegration 


Liper vals LoL 
turbulence 


(Ci Stee ee tee, i ile ee Giri la Ee Gee se eee ee Ei a ee se a fe cee Ste ls i, 


ee Oe ee ee ee BPP 


| 
| 
| 


SIT AND MITFJTI functicr alyorithmically the saue as 
fmeectmren Migk Gs Togeticr, tney compute the jitter WIF array. 
HiFJII is calied by JIT for each radial increment. The 


jitter MIF is given try 


pid 
M(B) = ona(* ee (3.29) 
- a: | 
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a 


wnere 29 iS a user ¢,eCitied input. AS with keaz guaiity, 


5 


this cuartity 1S appiied ac the aperture. 








de ee ee a 
| 
DAELE cae 
JII1 and MTFJIT Erogram Variables and Definiticrs 
{ 
Veliaris Portran Name See 

sup- Jain Derinitica 
i = 
| k (aj ' CK > 27/r | 
j Enis) E J 1 ses Janie beetle | 
i i 
| | 
Mn a st ee es ee mc ae ees oe emcee ne nee es i ne cr 

14, ERTIES 
This subroutire calculates tnat part of therial 
rklocring net dependert on the slant patn of the beam. The 
eguaticn used by GUTSAVG for the ;hase distortion due tc 
tnergal cleocming is 
Ad (x,y) = 27 Be n-1 "ry ax'] 
h 
tm 
: (2.20). 
Q a(h)exp[-secéJ/ a(h)+o(h) dh 
x h 
1S 
a oe ee 
ie Vocos€é + wh 


Note that tte first term in this expression is in variant 
with respect to keam sath while the second term contains 
such y-atk dependent variables as wind , wind due tc slew, 
and extirction ccefficients. The first term is evaluated cry 
this sukroutine and the remainder evaluated later in the 
frogranz inside the shcct angle iteration loop by subrcutize 
AV. 
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Standard valves used Icr eels “ay ana.’ are 
fee Ot teUiIx 10> Jyis, and 1.4, resoectively. 4 phase 
Pa emeere (leer isoeeGostrucited by Lteration or the irtersity 
intecral. Ncte that [frcgram performs this integration cver a 
half piace cf the aperture defined by -k<y<b and O<x<k where 
k is tké aperture diameter. The results are adjusted later 


Mrmestkrcttire PHVAR ficr this rethcd. 


JABLE XLIil 


PEIPHS Procram Variatles and Definitiozs 


Vériatle Fortran Name ee 
Due / Hain Definition 

ker) OK - 2m1/A 
Game t. - Refractive 


index tern 
(awe x 10g) 


Y const. zs it OmO beecie> cl Lic 
aeats (1.4f 


es Gone ts = Atmos. pressure 
(1.07x10S) Nyme 
GCin, ¥ ) Pai) PH{1,)) fnermal Lloouaing 


Ritts Milan Baten ie tlie Eee, ORs Sate RE AEE Pete, Meee ee HE Bann Satemer Betis Sel nan Sm i ee oe, 
oan 


(OnKS, —Dte ee e Oe 3 Me Bi I ee Saale tite Bae ime — Gicemie Glatt Sy ee pe Ce 


P 
| 
| 


18. RESFHS and CCEFFS 


Cnee the phase screen has been constructed by 
Ssubrcutine FRTPHS, sukroutine FESFHS removes the tear, tilt, 
and fccts curvature. Tke result is the poase aberraticn dve 
to thermal Elooming alone. Zernike polynomiais are utilized 
for e€xpressing these phase distortions with subroutine 
COEFFS froviding the required ccerficients. These pelync- 


Dialsé ale 
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i 
fo 


2, (a) ee (mean) (25 S30 


21 (X,Y) = aixX rx Cle) (3. oe 
2 {X,¥) = ay ioe oi oe) (3.2 
Z3 (X,V) = a3 (x*+y*) +a, (focus) (3. 203 


The €xfarsicn coefficients are given below and arc ccempeimes 
relative to a unifomr apenture aieighweng fimebion WY ine 
W(X, y) = 1 anside the aperture and zero elSewhkere. The iute- 
Gral¢ kave céen nultiplied Ly Zz to compensate for the half 


plane integration of ==uc rou Lincmemn no 


f,= <a, O¢(X,y)W (x,y) dxdy (3 Sp 
Ti= <4, o(X,y¥)W (x,y) xdudy ( doar Gy 
f,= 2a, o(X,¥)W (x,y) ydudy (3.238 
f= < (a ee (x,y) W(xX,y) dxdy (3. 38) 


the final p~kase correction is given by 
o(x,y)=o(x%,y) -£121 (x,y) -£222 (x,y) -£ 523 (x,y) 


(3.28 
ang no \ aay ) 


For a discussion of least Squares fitting ci Zeriike peilvaee 


tials, see jRef. 36]. 
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ioe eoscevartanco due £¢e tlermal plocning is 
computed Ly this subircutine using the residual phase screen 
crovided Ey RESPES. \Ncrralized with respect to the aperture 


field, tre variance 1s given ry 


SC” _ wane yim" (x,y) + U(x,y) o(x,y) (Gee 
UNG ny) 
ae aria a = ae 
TABLE XIY¥ 


PHVAH, HESPHS, ana CCzFS Program Variable Definiticrs 


| 
| { 
q 

| 
| 

{ { 
{ Variable Fortian Nae 
Sub” elaan Bet eo 
b(X, ¥) Pei ee Ht) Pnaase array j 
i Z . i 
SiG S20 SiGse0 Phase variance | 
. ; | 

A , Aj- A4 - Z2ernixe coetfrs. i 
{ 

| = oe PHMEAN = The exparsiocr } 
{ Pere eK oa coeificients j 
j eer Y es relative to a i 
j PHeCeuSs = UNLIZTOCM weignting \ 
7 ore seeye 
P Z 

Pies DV 


AY évaluates the path dependent part (second tern) 
ci tre trervai bloomirg eyuaticn 35.30 . (see Ssecticn FRIPES) 
GO(h) and O(£) are the absurgticn coefficient and scattering 
Seeiidelent, respectively, V¥, is the user specified wind and 
is tke slew rate. The denominator represents the tctal tran- 


sversée wind cOmpcnent across the ‘team. 6 is the zenith angle 
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and Sis the angle cf attack cf V, with respect to the beau. 
V has keen assumed tc be Opposite in direction ana paralieg 
to the target to trarsemitter LCS Slew motion. kote “aikeae 

that V , wail be appiied as a ccrstant the entire lencth of 


tne kéeam path. This is discussed in Chapter one. 


| 
| 
| 
| 


eS Ee eS 


TABLE AY 


AV Prograz Variables and Definiticns 


i 
| | 
ae | 
Variance Fortran Name { 
Sul ~ Hain vefinition 
| a{h)+#o(h) ALE - Total vextins tle i 

coefricient | 
l - | . | 
a (b) (ka) ” ALS = Anpsorpticn coerf. 
j séc (6) SECEG SECC FG secant or zenith i 

angle 
| V, (m/sec) VO VO Wind | 
Piri Ge PHILCCT Slew rete 
b (meters) HAT#C HATIMC feiyht of the j 
| atm atmOspkere | 
b (7é ters) oe HRA NG Helght of tre { 
j t transmitter 1 
|, ee a _-= 

lee  ELeGs 


The purpcese cf BLOOM is tc provide Treason atie 
results fer theraal Elccming dégradation when tke Irel 
values are btelow appicximately 0.3. In this region, the 
exponential Strehl relation predicts unacceptabiiy severe 
results. Therefore, sten the phase variance is greater than 
1.2, ELCCM computes an Irel value based on GUISMTF results. 
GUTSHIF is a full wave oftics propagation ccde using FFis 


ek 


Curve rat cclynopials were develczped using GUISHTF results 
Por Ecth a UHiLorm ara a truncated Gaussian aperture distri- 


fmerzen., Ihe Pesultire pPciyncmieis are 





Irel = * (3.41) 
“202005 + 2.914296 + .17230? 
il 
2 eceeninemmnaenccneetiea Piya 
Meee 2 64916 + 4.096030- (2. 42} 


ii the ;refile under consideraticn is a truncatelil Gaussiar, 
i.e. the aperture diareter greater than tne waist diameter 
cx tre kean, then the reutine uses tke truncated Gaussian 


polyncnial. Otherwise, a combination of the two are used, 
= e 2 2 (32743) 
Trel,, E (ays) |rren, + (d/b) eee 


auere d is the ayerttre diameter and b is the waist cCiazeéter 
Suttle Lean. [Ref. 27} 
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emma mm Samm mm es a ee me ee 


DR Ee ae 


EIOOM Erogram Variables and Definitions 


{ : 
{ l 
| | 
a | 
| Veriakie Portran Hane ce... | 
Sub Maan Definition | 
| o* Ss SiGee Phase variaace ! 
i . 
] Irely, at TELCCH Total blooming | 
| prcduced [rel { 
i { 
j Irely, TRU = Curve f£1t oles | 
{ LOE ia: Ort | 
| aperture dist. 
Itélig TRG = Curve fit Irel | 
Tor Gaussian 
| aperture dist. ! 
| | 
Os ce se a A ES NS SS SS AD ES SS ES SA SA AOD SD SD cD cD ee 
lS gees 


KIFAIM computes the atmospheric MIF and Lf Seceme 
fied, applies a tilt due to turbulence correction. ihe wee 
is giver Ly 


5 


= ; ] 
M.(D) = exp -(o/0) “3 | 1-a(5/p) 5] (oa | (2.44) 


where [ is the aferttre dianeter,o) is the coherence dian- 
eter and d is defined as 1-ALAE. ADAP is the fracticnual 
residuéeji tilt due to turbulence and is a user specified 
Farameter. If ADAP=0, M will represent the fully turrulence 
tilt ccrrected MIF. If AD@P=1, M will be the uncorrected 
MTF ard way ke writter as 


Sy, 
M, Co) = exp|-(B/00) | (3.45) 


ALAP way Le any value Eetween O and 1. 
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——— er ee ee Se ee ee ee ea ee ee Sas el ae ee eee come ee 


TABLE XVI 


MIFATM Procram Variatles and Definitioms 


Véeriacile Fortran Name 

— e sub Has Definitica 

B0) 1g fee Atmospheric ATF 
O9(MEtEIS) RHO 5 HO Conerence dia. 
CP (meters) D DIA Agere ui 6a24 . 


i ee 


| 

| 

| 

. 
2 —— | 
i 

( 

| 

| 

i 

| 


' 


> EP ~  EP i aP een e cie mmeS ee 


— 
CF 
—_ 


Zee LAE o 


t~ 





TREHEL Ccmputes the appreximacre Irel for a tilt 
corrected system. This value is used for ccmfarison with 
the Irel preduced by a prefect adaptive optics comjensatec 
system. The purpcse cf this ccrparison is to determine if 
noise will degrade tke AO system to an extent as tc make a 
full AC system undesirable. If this is the case, the [rcgran 


Meeiedy yay tilt correction only to the Lean. 


Subroutine FALRHO calculates the coherence diameter 
that wteuld result in the Icel value acaieved by a pertect AO 
system. This is done by successive calls to subroutine 
TRBREI and iterating p,- This new P, then becones a factor 
in the AC ccmpensaticr. Specifically it will cre used Ly 
Subrcetine MIFATM tc produce the atmospheric MTF of the AC 
corrected systen. 


84 


TEERFL and FSDRHC Program Variables and Defimiticns 


a eS E—————_———E_———————————— a 
| TABLES Via 

| 

| 


Variarie Fortran Name a... 
Sub azn Definition 
rp) Eee = AEMNOSPHE L1G ae 
Ma (?P) F1 IRRMIF Aperture MTF 
Trel ee EXD Irel orf perfect 
| AO systen 
a 2 





27. TEI 1S0 


The call to TITISO is wade when adaptive oftics nav 
not keen user specified ana iscrlanatic caiculaticns have 
Lot keen inribited. Also,if the signal to noise ratic at tk 
AO senscr is such that the AO system will provide tilt 
correcticn cnly, TLIIsSC is called , again, provided Zecyeeee 
Matic calculations héve not been inhibited. 

The purpese cf TLITISC is to include beam wancer dte 
to iscplarnatisgm. The 2-sigma-p tilt 1s Computed and 
combined with the jitter 2-siyma-p. This combined ter is 


then used tc compute the jitter MIF. 


dad 


23. FELUOT 





FELICT simply averages the reiative intensities due 
to a guitiplicative ard an RSS(roct sum Squared) approach t 
combining thermal blccming and the other propagaticn 
effects. 
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Le setae gg Mis MB ee es BO eet BH Se Be eis ee OE, 


© 


& 


O 


+ iL mee. 
Mee = ae = m (3.5416) 


2 


Tiel ais the result ci the -multiplicative approach 
m 


= = (Sd) 
Irel. Irel,, Trel. 


where Trel,,1is iemereieduew to (tice tare plocming and izel- is 


Mmmmeinre! due Ec Leam cuality, jitter, turculence , iscplara- 


Misi dgiG acaptive optics effects. Ihe 255 apgroach is given 


ry 
rel, .. “h(a eee (res! = (248) 
amo re" th 


| 
| 
| 
| 


TABLE XIX 


ee ee | 


RELTOT Precgiram Variatles and Definitions 


| 
| | 
i Vériatie Fortran Name 
| Sub ' Main Dorn On | 
ere Le a pO en. Totai Irel 7 

| 
| irel TR hoe Pee beh S | 
| Ss neue. 
| cel at Se gh O Ea Irel ae { 
aes icati ve | 
metho 
x L EXD Irel of perfect | 
AO systen 
| Irel, = HR OE eel die to” j 
] azl otner etfects 
j j 
Irel - MeL eG | Irel due to j 
| tb thermal blooming | 
| 

L 


i 
| 
| 
| 
| 
| 
is 
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24. ISOPLA 





This coutine calculates the MTF that characterizes 
the isoplanatic erfect on the predictive or "look-ahead" 
adaptive optics system. The look-ahead angle is the major 
input parameter to this routine. The Fortran code was 
written by D.L. Fried [Ref. 38]. Fried developes the 
isoplanatic dapendency of the AO system in terms of the 
effective antenna gain of the laser transmitter. The MTF 


that Fried formulates is given by 


TT 


—2-—5 
ae! =f rsa.san3% -sec(8)-u(vsec(6/p),¢)dod (3-49) 


where 


x =eyseca0 gay (3.50) 
and 


oP 


U(x, >) - fos faeom -§;[1+2(xh) cos + (xh) *| 
(3.51) 


5 
-} [1-2 (xh) cos + (seh) + | ‘tah 


6 is the zenith angle and y is the target lead angle. 


For detailed jdiscussion of the theory and explicit develop- 


ment of tue Fortran code, see [Ref. 39]. 
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HUMGCOR TLCS tLe 2er@gecra@er Bessel function based cn 
Pies i0EUt arguments this routane is Called py FARFLE it the 


@abreculaticn cf the far-£ieid irradiance. 
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NkrEIs (6, 180 (eee Zoe i), 1. is 1, 50) je ew lDid 
WRITE (6, 180 (OE Seno 1), C1 = 1, S0}e aA 
WRITE (6,185) (DESCE (27,11) ,11=1,50), N1 
Sorter (6, 185 DESGR (Zo 0b), L l=, 90), NZ 
mRiae (6,185 PE CGR oe, t= oo). NS 
WRITE (6,185 (DEECR 30,11) ,11=1,50) , N4 
Mike (6, 185 DESG Mia ly), L1= 1, 5035, NS 
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eee assess HATA REALIGNMENT SECICER == 2322222225 S2252S25 2525 2252>==2=€ 
c 
Meena. 2-0.) SCALSO-[1a/s. 
Meret WA VEEC.NE. G2) VARGD= (2. *PI*nAVEBC) **2 
SGJITO=SIGJIT 
LF Tse Eos i). AT AF=0. 
Mee CNOTSC.5EQ.1) ITSGAKO=1.E10 
coe so stad s 
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b pETMx= Peeina/ £1) ; 
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RAEVGD=HTSANS~-EGEND 
: IF (HABYGD.LT.~C.~) HAEVGD=0. 
R=REARTH+ETRANS 
RS=REARTH+HSAT 
NFLGSN=0 
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ee ee ee ee ee ee ee ee = ee —_—— <r ee ee ee ee ee ee ee ee ee ee es ee ee ee ee 
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‘S ee ee ee ee ee ee ee ee eee eee ee ee ee eee ee ee ee ee ee ee oe 
Cc 
ANGOFP=LCFF/SR& 
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& 
eee ne nn a nn sa ss se nn nnn ee ee ee een ee ee 
c m@neuleE COORDINATE IBANSLATIONS DUE TO CFF=SET 
Cc — ee Oe Oe ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee es ee ee we ee ee ae 
iS 
LY=+2. *R*SI Moss) ee 
LZ=2.* R*¥ SIN 
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c a me ee ee ee ee ee ee ee ee ee ee ee oe oe oe a eee me ee ee ee ee ee ee ee ee ee ee ee ee ee oe 
S COMPUTE EARTH CENTER ANGULAR aAATE 
ec ee ee ee ee ee ee ee ee ee ee ee ee ee ee eee ee es ee ee ee ee ee ee ee ee ee ee ee ee ee 
ec 
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Cc — a oe oe eee ae =e a oe oe ee oe oe oe ee ee ee ee ee ee ee ee ee ee ee oe eo ee ee ee oe ee ee i ee ee ae eae 
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Cen rr ne ne nn we ee ae een ee ee ne ee ee een ee — a oe eee eee ee ee ee ee we ee eee 
€ MeaU5St IRTENTANECUS SLEWRATE FCR ROTATICN OF EARTH, 
Cc assuming broth are cclinear. 
——— ee ee ee ee ee ee es ee ee 
e 
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oo — 2 ae ee ee ee ee ee es ee es ee ee ee ee en es es ee eae ae a ae ee ee 
c IF SATELITE ALZTILIUDE IS 48 TEE OCOkDE SRO eee soe 
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C INTEGRATICN BEICW. HENCE WE CEFINE DX NCh. 
Cc 
DX=DIA/NA 
C 
EC== ce ee ee ee ee we ce ee ee ee ee ec ee ee ee ee ee eee eee ee 
C THE CALL TO UOCST DETERMINES TEE CONSTANT THAT MAKES TRE EXIT 
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Git Doce” (PICTAL, Bi) 
CALL PARFID (D1,115C,IRRMTF,ML,NY,HSATD, WAVE, PTOTAL, PMAXO) 
G 
C= — <<a ee we ce ce ee ee ee es ew ee ce ee es se es ee ee es ee es ee < e s e es e 
E IF THE WAVEBQ IS NCT SPECIFIEL AS INPUT, COMPUTE ON BASIS OF 
c TDFLMT. 
ce ee ste i 
IF (WAVEEC.EQ.C.) CALL DETWAV (IRRMTF,DX,N&4) 
C TF IDPLAT EAS Sel SEECIFI lO), CCreUor mye 
Cc 
IF {TDPLET.NE=C-) GO 1 
CALL BQIEEL (IFEZTF, 4,DX,7BQ) 
TDFLST=SCET(1-/71EC) 
jo CONTINUE 
Cc — oe ow cee OP ee ee ee ee ee ae ee es ee ee ee es es ee es es ee ee es ee ee ee a es es es i es es es es ee ee ae i 
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C 
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Ce ce ce ee re ce ce ccc cc ce ec ree ce ce ee ee ce ee ee ee eee ee ee ee ee ee 
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COMPUTES THE ITEERMAL BLOOMING EHASE VAERIANCE 
VARIANCE IS CCHFOTED LIKE THE STREHEO RATIO SeCOREING Tomas 
FIFLD AS A WEIGHTING FUNCTION. 
COMMON /PEBFLL/ II}, CIAOES, BEAMSZ,U0 
DIMENSION PH(#LE,#LS 

RCSQ=(DIAayz.) #92 

RISQ=(DIACSS/2.) *4*2 

MLM=ML-1 

MIHM=MLH-1 

DELN=DIA/PLM 

XI=DELN/<. 

YI=-DIA/2.+XI 

X=XI 

SUM1=0. 

SUE2=0. 

SUM3=0. 

DC 30 I=1,MLHM 

y=yI 

XSC=x*X 

DO-20 J=1,MLé 

BSC=XSOQ+Y4Y 

Tre (hS0.G7.8CSC. OhanSCeltoRT Sc mccm Cmte 
CALL PIELL (SCT (ESC) ,0R) 
SUM1=SUM1+0 

SURZ= SUN 4 +DRSER (12d ** 2 
SUM3=SUM2+UR*EFE (Id 

CORTINUE 

Y=Y4+DELN 

CCNTINUE 

X=X+DELN 

CONTINUE 

-- ADJUST FOR HALF ELANE INTEGRATION 
SUMI=SUM 1+2. *LEIN**2 

SUM2=SUMZ+2. *DEIN**2 

SUM3=SUMi+2. *DELK**2Z 
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SiGSQO0=SUM2Z7 suet 5sUN3/7/SU8 1) **2 
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SUBROUTINE TREREL (IRPMTF,DIA,N4,DX, TREL, EHO) 
B---— THIS SUBHCOTIN? CCMECTES THE RELATIVE INTENSITY DUE IG A TOTAI 
C---— TILT ORRRECTEL SYSTEM FOR TUREULENCE CNIY. 

BEAL ae 

PI=2.*ARSIN(1. 

SUM=0. 

DX 2z=Dxi/2- 

X=DX2 

Demo 1=1,N4 

FASIREMTE {7) 

@kiL MIFAIM (%,CIA, RHO, 0. ,F2) 

SUN=SUM+F1*Fz2* 

X=X+DX 
10 CONTINUE 

TREL=SOM*2.*E19DX 

RETURN 

END 
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Cm----THIS SUBROUTINE CCMNETIES THE BESITDUATL TILT ous 2G sor tani ise 
c 
COMMON /ATMO/ BATHO 
N=100 
F56=5. /6.- 
Eee nCn: 
DS=(HATMC-HTRABS) /N 
H=BABVGD+IS/2 
SUB=0. 
G 
DO 10 I=1,N 
C 
- CNLL (CNZE (nh Che) 
: Z=SECOMG*H 
ARG1= D+Z*OMX a541 (oooRe **2)** E56 
ARG2= D-Z*OMX Veet LORY 22] 28236 
ARG3=(2*4*2* (OM Y**Z+OMY* #2)) **# FSG 
ARGY=D** E53 2 
Tpee ee eons S45 hee egH Teo Reese 
ARGoO= Z*CMX) **24 (D*#Z*OMY) **2) **# E56 
cj 
Bee a eee ee eee eee eee al eae 
=n? 
ue CONTINUE 
eaara THIS IS TBE REEULIANIT 2—SIGKA=EF TILT DUE TO 15 0clLANaALLo 
P RESTLT =SUM*DS*SECCNG*2.91/D** 292. 
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